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Abstract 

Theoretical studies of the tunnehing current and emission spectrum of a single electron transistor (SET) 



tinder optical pumping are presented. The calculation is performed via Keldysh Green's function method 
Q within the Anderson model with two energy levels. It is found that holes in the quantum dot (QD) created by 
'""'optical pumping lead to new channels for the electron tunnelling from emitter to collector. As a consequence, 
an electron can tunnel through the QD via additional channels, characterized by the exciton, trion and 

o , 

y—^ oiexciton states. It is found that the tunnelling current as a function of the gate voltage displays a series of 

o 

\^ sharp peaks and the spacing between these peaks can be used to determine the exciton binding energy as 

I 

> associated with the electron-hole recombination in the exciton complexes formed in the QD can be controlled 

a 

^ both electrically and optically. 

o 
o 
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I. INTRODUCTION 



The tunnelling current characteristics of a single electron transistor (SET) have been extensively 
studied and interesting physical phenomena, such as Coulomb blockade and Kondo effect have been 
found. The main structure of an SET consists of three electrodes (source, drain and gate) and 
the active region consists of one or more quantum dots (QDs). Although the physics of Kondo 
effect is more profound than that of Coulomb blockade, the main device application of SET is based 
on the Coulomb blockade effect. The tunnelling current of an SET is sensitive to the size of the 
QD. The charging energy, AC/, is typically larger than the energy level spacing, A£', caused by 
the quantum confinement in large QDs. Consequently, the tunnelling current displays a periodic 
Coulomb oscillation with respect to the gate voltage. This feature of periodic Coulomb oscillation is 
a consequence of homogeneity of AU and negligible AE. However, such a Coulomb oscillation was 
observed only at very low temperature due to the small charging energy. Note that the charging 
energy must be greater than the thermal energy and tunnelling rate multiplied by h in order for 
an SET to operate properly. For apphcation in high temperature SETs, nanoscale QDs with large 
charging energy are desired. 

To date two major kinds of nanoscale QDs have been used to construct SETs: the Si/Ge QDs and 
the InAs/GaAs self-assembled quantum dots (SAQDs). Silicon (Si) or germanium (Ge) QDs can be 
miniaturized to reach nanometer scale with advanced fabrication technology. Nanoscale Si/Ge SETs 
typically display non-periodic Coulomb oscillations at room temperatures, as a result of the large 
charging energy (AC/) and the unequally spaced values of AE, which are comparable to AU. Even 
though Si/Ge SETs can operate at room temperatures, it is difficult to gian full understanding of 
their tunnelling current characteristics because of the multi-valleyed nature of Si/Ge and the unknown 
Si/Ge-Si02 interface properties. Furthermore, the electron-hole recombination rate in Si/Ge QDs is 
not large enough for useful application in optoelectronics, since Si and Ge are indirect-band-gap 
semiconductors.^^ 

The InAs/Gas SAQDs, on the other hand, are better understood, and because InAs is a direct- 
band-gap semiconductors, SETs constructed from InAs/Gas SAQDs can have useful apphcation 
in optoelectronic devices such as single-photon generators and single-photon detectors. Recently, 
InAs/GaAs self- assembled quantum dot (SAQD) embedded in a p-n junction has been employed to 
generate single-photon emission for application in quantum cryptography. Another way to achieve 
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single-photon generation is to use isolated SAQD with optical pumping. Such studies have received a 
great deal of attention recently. The spontaneous emission spectrum of a QD typically exhibits 
coexisting sharp emission peaks, which have been attributed to the electron-hole recombination in 
the exciton, trion and biexciton formed in the QD. The antibunching feature of the emitted single 
photon in this kind of device has been demonstrated. ^^'^"^ In addition, gate electrodes have been used 
in some experiments to manipulate the emission spectrum.^^"^^. 

Many theoretical efforts based on the Anderson model and Keldysh Green's function method have 
been devoted to the studies of transport properties of SETs, including the dc and ac tunneUing 
current. The studies on the ac tunnelling current of SET mainly focus on the photon-side band 
phenomena caused by the microwave pumping"'''^. Only a handful theoretical studies considered the 
optical properties of SETs. Recently, we have studied the intraband transitions of SETs made of 
InAs/GaAs SAQDs for application as an infrared photodetector with operating wavelength near 
10//m.^^'^^ On the other hand, the interband transitions in a QD have been adopted frequently to re- 
alize single-photon generation. Theoretical studies of the electronically driven single-photon generator 
(SPG) have been reported, ^''''^^ whereas theoretical studies on optically pumped SETs (including the 
coupling with leads) are still lacking. It is the purpose of this paper to provide a theoretical analysis 
for such a system. Furthermore, the measurement of exciton binding energy in a QD remains illusive, 
since in a typical photoluminescence measurement or the emission spectrum of a single QD, one can 
only observe the exciton recombination, but not the free electron-hole recombination. The exciton 
binding energy is an important physical quantity as it provides information about the effective dielec- 
tric constant and the charge distribution of electron and hole in the QD. In our current theoretical 
study, we will show that via the measurement of the tunnelling current of an SET made of a nanoscalc 
QD under optical pumping one can determine the exciton binding energy of a QD unambiguously. 

The schematic diagram for the SET system considered is shown in Fig. 1. Here a nanoscale 
QD is embedded in a n-i-n junction. This is different from the electrically driven SPG, where a 
QD is embedded in a p-n junction. Therefore, the physical process considered here is also quite 
diffrent. Using optical pumping to create electron-hole pairs in the SET, one can not only electrically 
manipulate the single- photon emission spectrum arising from exciton complexes in the QD, but also 
optically control the tunnelling current of the SET. In this paper, we will consider both the hole- 
assisted tunnelling current due to optical pumping and the electrode-controlled spontaneous emission 
spectrum. 
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Our calculation is based on the Keldysh Green's function approach^^ within the Anderson model 
for a two-level system. A effective-mass model, which takes into account the anisotropy and inhomo- 
geneity, is used to estimate the inter-particle Coulomb energies. We find that the optical excitation 
creates holes in the QD, which provide new channels (via the electron- hole interaction) for the electron 
to tunnel from the emitter to the collector. As a consequence, an electron can tunnel through the QD 
via four additional channels, characterized by the exciton, positive trion, negative trion, and biexci- 
ton states. Each addition channel can generate a new plateau (or oscillatory peaks) in the tunnelling 
current characteristics in addition to the typical plateau (or peaks) caused by the Coulomb blockade 
effect. This gives rise to a rich tunnelling current characteristics and it can be used to determine the 
exciton binding energy. 

This paper is organized as the following. In section 11, we derive the tunnelling current of SET under 
optical pumping. In section III, we calculate the polarization needed for describing the spontaneous 
emission spectrum. In section IV, we present calculations of the inter-particle Coulomb interactions 
within a simple but realistic effective-mass model. In section V, we discuss the results of our numerical 
calculations on the tunnelling current and spontaneous emission spectrum and demonstrate that 
the exciton binding energy can be extracted via the measurement of tunnelling current. Finally, a 
summary and concluding remarks are presented in section VI. 

II. TUNNELLING CURRENT 

The system of interest is shown in Fig. 1, which consists of a single InAs/GaAs self- assembled 
quantum dot (SAQD) sandwiched between two n-doped GaAs leads. Electrons are allowed to tunnel 
from the left lead (emitter) to the right lead (collector) under the influence of an optical fleld. The 
Hamiltonian for the system is given by 

a,i=e,h ]i,a;£=L,R 

£=L,R;k,a 

+ E(Aoe^'^°*6e,k,.Vk,-. + >^le-^''%K~A,Ka) 

k,fr 

i,k,q,CT 
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Q 

where the first term describes electrons in the QD. We assume that the quantum confinement effect 
is strong for the small InAs QD considered here. Therefore, the energy spacings between the ground 
state and the first excited state for electrons and holes, AE^ and AEh, are much larger than the 
thermal energy, ksT, where ks and T denote the Boltzmann constant and temperature. Only the 
ground state levels for electrons and holes, and E^, are considered in the first term. The second 
term describes the kinetic energies of free carriers in the electrodes. Note that in the current setup, 
the gate electrode does not provide any carriers, but merely controls the energy levels of the QD. The 
third term describes the coupling between the QD and the leads. Note that only the conduction band 
of QD is coupled with the electrodes. The fourth term describes the interband optical pumping with a 
frequency Uq, which is in resonance with the transition energy for an electron-hole pair in the wetting 
layer. We treat the electromagnetic field as a semiclassical field. The fifth term describes the optical 
phonon assisted process for electrons in the wetting layer to relax into the InAs QD. The last term 
describe the coupling of the QD with the electromagnetic field of frequency cu. X= —/ir^ is the Rabi 
frequency, where //^ =< > is the matrix element for the optical transition and £ is electric field 
per photon. a^^(A^) and a (A) denote the creation and annihilation operators of a photon (phonon), 
respectively. 

Due to the large strain-induced splitting between the heavy-hole and light-hole bands for typical 
QDs, we only have to consider the heavy hole band (with — ±3/2) and ignore its coupling with 
hght-hole band caused by the QD potential. Thus, we can treat the heavy hole as a spin-1/2 particle 
with a —l,^ representing = ±3/2. This treatment is convenient for algebraic manipulations in the 
calculation of Green's functions. Because the effect of inter-particle Coulomb interactions is significant 
in small semiconductor QDs, we take into account the electron Coulomb interactions and electron-hole 
Coulomb interactions by adding the following terms 

Hj = Ui^idl^^di^crdl^_^di^-cr + J3 Uijdl^^di^fjd^j^_^dj^-„, (2) 

i,a i^j;a,a' 

where dj^ and di^^ are creation and annihilation operators for particles with spin cr in the ith energy 
level of the QD. Once the Hamiltonian is constructed, the tunnelling current of SET can be calculated 
via the Keldysh Green's function method. We obtain the tunneUing current through a single dot 
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(see appendix A) 

where fii^) and /^(e) are the Fermi distribution function for the source and drain electrodes, respec- 
tively. The chemical potential difference between these two electrodes is related to the applied bias 
via IjLl — piR — eVa- and Tii{e) denote the tunnelling rates from the QD to the left (source) and 

right (drain) electrodes, respectively. For simplicity, these tunnelling rates will be assumed energy 
and bias-independent. Therefore, the calculation of tunnelling current is entirely determined by the 
spectral function A = ImGl ^.(e), which is the imaginary part of the retarded Green's function ^{e). 
Eq. (3) is obtained provided that the condition ^l{r) >3> 7e,c and ^l{r) ^ Reh are satisfied, where 
7e,c is the captured rate of electrons. Because of the phonon bottleneck effect, we expect the rate for 
the electron to relax from the wetting layer to the QD is small. Thus, we have ignored the terms 
involving phonons in the calculation of the tunnelling current. Besides, photon current has also been 
neglected, since the electron-hole recombination rate, R^h is much smaller than the tunnelling rates. 

The retarded Green's function, Gl^{e), can be obtained by the equation of motion of Gl^{t) = 
— i6{t) {{demerit), dl^{0)}), where 6{t) is a step function. The curly brackets represent the anti- 
commutator and the bracket (...) denotes the thermal average. The Fourier transform of Gl^^{t) 
is given by 

g:A') = r dtGiAty^""''^' (4) 

J —oo 

with 1] being a positive infinitesimal number. G^ ^ is the full Green's function, which includes all type 
of interactions. Given the conditions ^l{R) ^ 7e,c and I^LiR) ^ Reh, we can drop the fifth and the last 
term of the Hamiltonian in the calculation of Gl^^{e). To solve G'g p.(e), we consider only the lowest 
order coupling between the electrodes and the QD. The equation of motion for Gl^^{t) leads to 

{e-E, + z^)G;,(e) = 1 + U^Gl^ie) - U,h{G\A^) + G:^,{e)), (5) 

where the two particle Green's functions, Ggg(e), Ge^i(e) and Ge/i,2(^) arise from the particle corre- 
lation and they satisfy 

(e - {E, + Ue) + ^^)GL(e) = N^,-, - U,n{Gl^,{e) + Gl^^{e)l (6) 
(e - {E, - Ueh) + iy )G^w(e) = + UeGl^^^e) - U^uGluhi^), (7) 
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and 

(e - {E, - Ueh) + iY)G:h2ie) = nh,-a + UeGl^,^{e) - U^hGl^hi^). (8) 

Ne-a, nh^a ^nd rifi-a in Eqs. (6)-(8) denote the steady-state electron and hole occupation numbers. 
The two-particle Green's functions are coupled with the three-particle Green's functions defined as 

Glheli^) = {Ne,-anh,ade,a,dla.), Glhe2{^) = {Ne-anh-ade,a,dla) and G'eM(e) = {nh-anh,ade,a, dl^.) . 

The equations of motion of the three-particle Green's functions will lead to coupling with the four- 
particle Green's functions (two electrons and two holes), where the hierarchy terminates. Thus, these 
three-particle Green's functions can be expressed in the following closed form 

+ !^h2zi ) (Q) 

e-{E, + U,-2U,n)+i\'' ^ ' 



1 - nh,a 



e-{Ee + Ue- 2Ueh) + i 

and 



rr)' (10) 



1 _ jV _ 

H ^SzZ ) (II) 

Eqs. (9) and (10) describe the mixed amplitudes for the propagation of an electron either in the 
presence of another electron (with opposite spin) plus one hole, or another electron plus two holes. 
Eq. (11) describes the mixed amplitudes for the propagation of an electron either in the presence of 
two holes, or two holes plus another electron. Substituting Eqs. (9)-(ll) into Eqs. (7) and (8), we 
obtain, after some algebras, the retarded Green's function of Eq. (3) 

^e,a\^) - U - ^^e,-a)\ ^TT^IV 

e — Ee + 

e- Ee + Ueh + i'^ e-Ee + 2Ueh + 

j,r r 1 ~ ^^h,a + nh-a) + nh,anh-a 

^ ^'-^^ e-E^-U^ + i^-f 
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e- Ee-U^ + Ueh + i^ 

I nh^ijUh^-cr I 

e-E,-Ue + 2Ueh + i^ 

Here Gl^{e) contains an admixture of six possible configurations in which a given electron can 
propagate. These configurations are: empty state, one-hole state, two-hole state, one-electron state, 
one-electron plus one-hole state, and one-electron plus two-hole state. In Eq. (12), Fg is the electron 
tunneUing rate Fe = F^, -|- F^^. The particle correlation has not been included in the tunnehng 
rates of Eq. (12). Such approximation has been known to be adequate for describing the Coulomb 
blockade effect, but not the Kondo effect. The electron occupation number of the QD can be solved 
self-consistently via the relation 



Ne. 



f deTLfiie) + TRfR{e) 

- = -J V fTTr^ ^'^^ 

A^e,(7 is hmited to the region < A^e,^ < 1- Eq. (13) indicates that the electron occupation numbers of 
the QD, Ne-a and Ne^a, ai'c primarily determined by the tunnelling process. To obtain the electron 
and hole occupation numbers {rie-a = ne,a and Uh-a = nh,a) arising from the optical pumping, we 
solve the rate equations (see appendix B) and obtain 

7e,cA^e,k(l - K) 



and 



(14) 



";/),c'^ /),k |.N 

l^h — ^h-(T — ^h,(T — , -D / , AT \ I ^ \ ) 

lh,cNh,k + Rehirie + Ne) + Fft,^ 

where 7e(/i),c and A^e(/i),k denote the captured rate for electrons (holes) from the wetting layer to the QD 
and the occupation number of electrons (holes) in the wetting layer. See Eqs. (37) and (38), Nf.(h),k 
is a linear function of Pexc for low pumping intensity and small 7e(/i),c- T/i.s denotes the non-radiative 
recombination rate for holes in the QD. Because 7e,c^e,k/re 1, rig has been neglected in Eq. (12). 



III. SPONTANEOUS EMISSION SPECTRUM 



When electrons and holes appear in the QD, their recombination leads to single-photon emission. 
The optical polarization of the emitted photon depends on the spin polarization of the electrons. 
An electron with spin +1/2 (-1/2) can recombine with a heavy hole of angular momentum -1-3/2 
(-3/2) to emit a circularly (cr+ora^) polarized photon. In some experiments, up to four coexisting 
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peaks (associated with exciton, negative trion, positive trion and biexciton) have been observed in 
the emission spectrum. The emission process is described by the last term in the Hamiltonian 
given in Eq. (1). The bias-dependent emission spectrum can be obtained by finding the nondiagonal 
lesser Green's function at nonequal time/ To simplify the problem, we assume that the applied bias 
is large enough to create the biexciton state (two electrons and two holes) in the QD. 

With the above assumption, polarization P(a;) = G^fJ^uj) can be calculated by the equal-time 
Keldysh Green function, Gff^{t,t) — —i{a'^{t)dh-a{t)de,a{t)) (or density matrix method). Prior to 
the calculation of the polarization P{u!), a unitary transformation has been used to remove the 
phase of the optical transition term. The renormalized energy levels of the electron and hole become 
ee = -Eg — f and = Eh — ^. Furthermore, the hopping terms between the leads and dot, V^=L,_R;k(^) = 
Vj=L,R;]iexp^'^ 2 ' become time dependent, in which the energy and time dependence of the coupling 
are factorized. This factorization leads to time-independent tunnelling rates and it simplifies the 
calculation. Solving the equation of motion for Gff^{t,t) — —i{a^{t)dh,-cr{t)de,a{t)) and considering 
the steady-state condition, we obtain the lesser Green's function 

^< ^ -A.F - (Ue - Uek)Vi{u;) - (Uu - UeH)V2(u;) ^ ^^g^ 

Q + Cft ~ Ueh + 

where 

= ($|ata|$)(l - 7Ve,a - nh,-a) - N,^anh,-a, (17) 

r = Fg + F/i, which leads to a broadening of the spectrum, and |$) denotes a photon state. Vi{u!) 
and 7^2 (<^) are given by 



7^i(a;) = -AWe,-.( ^ """" 



+ eh + Ue- 2Ueh + ?f 



+ p) (18) 

^e + ^h + Ue + Uh- SUeh + h 



and 



1 - N.-a 



ee + eh + Uh - 2Ueh + 



ee + eh + Ue + Uh- 3Ueh + 



(19) 
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The first term in the filhng factor T [Eq. (17)] arises from the stimulated process, which vanishes in the 
vacuum state. The second term, N^^unh-c is due to the spontaneous emission process. This is a quan- 
tum effect of the electromagnetic field. Expressions of Eqs. (18) and (19) are derived, respectively, from 
the equation of motions for —iia"^ {t)Ne-a{t)dh,-a{t)de^a{t)) and —i{a){t)nh^„{t)dh-cr{i)de,a{i))^ which 
are coupled with -i{a\t)nh,a{t)Ne-a{t)dh-a{t)de,„{t)). The term -i{a^{t)Ne-a{t)de,a{t)dh-a{t)) 
describes an electron-screened electron-hole recombination process in a negative trion. Similarly, 
—i{a^ {t)nh^a{'t)dh,-a{'t)de,a{'t)) , describes a hole-screened electron-hole recombination process in a posi- 
tive trion, and —i{a^{t)nh,a{t)Ne-a{t)de,a{t)dh-a{t)) describes the biexciton-to-exciton transition pro- 
cess. In the photon vacuum state, we obtain 



V{u)/{N,^,nh,-.) (20) 



^ Al - N^-^){1 - nh,a) , {l-nh,a)Ne-a 



+ 



Eg-Ueh-oo + iV/2 Eg + Ue- 2Ueh - uj + iT/2 

(1 - Ne-a)nh,cr 



Eg + Uh- 2Ueh -0J + iT/2 

^ Eg + Uh + Ue- 3Ueh + iT/2^- 
The spectrum of the polarization V{u)) displays four peaks, which are attributed to the exciton 
X, negative trion X~, positive trion X+ and biexciton (biexciton decaying to exciton). The 
corresponding peak positions occur at u = Eg — Ugh, = Eg + Ue — 2Ueh, ^ = Eg + Uh — 2Ueh and 
CO — Eg + Ue + Uh — 'iUeh, which are significantly influenced by the particle Coulomb interactions. 

IV. ENERGY LEVELS AND INTER-PARTICLE INTERACTIONS 

According to Eqs. (12) and (20), particle Coulomb interactions will significantly affect the tun- 
nelling current of SET and the emission spectrum of single photons. To illustrate this effect, wc apply 
our theory to a pyramid-shaped InAs/GaAs QD. First, we calculate the inter-particle Coulomb inter- 
actions using a simple but realistic effective-mass model. The electronic structures of InAs/GaAs QDs 
have been extensively studied by several groups with different methods. The electron (hole)in 
the QD is described by the equation 
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= £;e(ft)'0e(?j)(r), (21) 

where m*^{p,z) (a scalar) denotes the position-dependent electron effective mass, which has m*^ = 
0.067me for GaAs and m*j = 0.04me for the strained InAs in the QD. mjl^(p, z) denotes the position- 
dependent effective mass tensor for the hole. It is a fairly good approximation to describe rri^{p, z) 
in InAs/GaAs QD as a diagonal tensor with the x and y components given by m*~^ = (71 -|- ^2)/i^e 
and the z component given by m*i^^ = (71 — 2^2) /"rrie- 7i and 72 arc the Luttinger parameters. 
Their values for InAs and GaAs are taken from Ref. 34. Vqj^{p, z) {yQ]j{p^ z)) is approximated by a 
constant potential in the InAs region with value determined by the conduction-band (valence-band) 
offset and the deformation potential shift caused by the biaxial strain in the QD. These values have 
been determined by comparison with results obtained from a microscopic model calculation^^ and we 
have Vq^ = — O.SeF and Vqj^ — — 0.32ey. The eFz term in Eq.(21) arises from the apphed voltage, 
where F denotes the strength of the electric field. Using the eigenf unctions of Eq. (21), we calculate 
the inter-particle Coulomb interactions via 

Ui,, = drj dv2 ^ V 22 
J J eo|ri-r2| 

where — e,h. nj(ri) denotes the charge density, eo is the static dielectric constant of InAs. We 
have ignored the image force arising from the small difference of dielectric constant between InAs and 
GaAs. Although the Coulomb energies are different in different exciton complexes, their difference 
is rather small.^^'^^ Therefore, only one set of Coulomb interaction parameters has been used in this 
study. 

For the purpose of constructing the approximate wave functions, we place the system in a large 
confining cubic box with length L. Here we adopt L — AQnm. The wave functions are expanded in a 
set of basis functions, which are chosen as sine waves 

tknlmiP, (t>, z) = ■^^sm{kix)sm{kmy)sm{knz), (23) 



where kn = mT/L,km = m-K/L^ki = in/L. n, m and £ arc positive integers. The expression of the 
matrix elements of the Hamiltonian of Eq. (21) can be readily obtained. In our calculation n — 20, 
m — 10, and ^ = 10 are used in solving Eq. (21). Fig. 2 shows the lowest three energy levels of the 
pyramidal InAs/GaAs QD as a function of QD size. The ratio of height and base length is fixed at 
h/b = 1/4, while h varies from 2.5 nm to 6.5nm. Diagram (a) and (b) denote, respectively, the energy 
levels for electrons and holes. We can see that the energy level spacing between the ground state and 
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the first excited state is much larger than ksT — 2meV, which will be considered throughout this 
article. Therefore, the lowest energy levels, Eg and Eh, are adopted in the Hamiltonian of Eq. (1). 
The intralevel Coulomb interactions Ug and Uh and interlevel Coulomb interaction Ugh are calculated 
using Eq. (22). Fig. 3 shows the intcr-particle interactions as functions of QD size. The strengths of 
Coulomb interactions are inversely proportional to the QD size. However as the QD size decreases 
below a threshold value (around b — 12nm), 17^ is significantly reduced due to the leak out of electron 
density for small QDs. These Coulomb interactions approach approximately the same value in the 
large QD limit. This indicates similar degree of locahzation for electron and hole in large QDs. We 
also note that Ueh is smaller than Ue in the large QD. This is due to the fact that in large QDs the 
degree of localization for the hole becomes similar to that for electron, while the anisotropic nature 
of hole wave function reduces Ueh- The repulsive Coulomb interactions, Ue and Uh, are the origin of 
Coulomb blockade for electrons and holes, respectively. The attractive Coulomb interaction Ueh gives 
rise to the binding of the exciton. To study the behavior of tunnelling current and bias-dependent 
spontaneous emission spectrum, we consider a particular pyramidal InAs/GaAs QD with base length 
b — 13nm and height h = 3.5nm. The other relevant parameters for this QD are E^^ — —O.lAeV, 
Eh,o = -0.125eV,Ue = IG.lmeV, Ueh = lQ.7meV and Uh = IS.bmeV. 

V. RESULTS 

In this section, we discuss our numerical results for the tunnelling current and spontaneous emission 
spectrum. For simplicity, we assume that the tunnelling rate — Fr — O.BmeV is bias-independent. 
We apply a bias voltage T4 across the source-drain and Vg across the gate-drain. With the apphed 
voltages, the QD electron and hole energy levels will be shifted to Ee + aeVa — fieVg and Eh -\- aeVa — 
jSeVg, respectively. In our calculation, we assume a = 0.5 and (3 = 0.7. Meanwhile the chemical 
potentials of the electrodes corresponding to a Fermi energy Ep = GOmeV (relative to the conduction 
band minimum in the leads) are assumed to be TOmeV below the energy level of Ee at zero bias. 
Other parameters adopted are: = 0.2meV and Reh — lO/ieV. 

A. Effects of optical pumping on the tunnelling current 

Applying Eqs. (3), (13) and (15), we solve for the electron occupation number Ne = Ne^a = A^e.-o- 
and the tunnelling current J. Fig. 4 shows the calculated results for Ne and J as functions of 
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gate voltage with and without photo-excitation at zero temperature and = 2mV. Sohd hne 
and dashed hne correspond to / = (no pump) and / = 0.9 (with pump), respectively. Here / 
is a dimensionless quantity, defined as / = 7ft,c^?i,k/r?i,s which is proportional to the pump power. 
The electron occupation number displays several plateaus, while the tunnelling current displays an 
oscillatory behavior. Both are typical behaviors due to the Coulomb blockade. The manipulation 
of gate voltage can tune the energy levels of the QD. For the no-pump case (7 = 0), carriers in the 
emitter electrode are allowed to tunnel into the QD as the gate voltage exceeds the threshold voltage, 
Vgs{Ri lOOmV). When the energy level of Ee + aVa — /3Vg is below the emitter chemical potential, 
but Eg + aVa — pVg + Ue above the emitter chemical potential, A^e (for a given spin) reaches 0.5 
and displays a plateau which is caused by the Coulomb blockade associated with Ue (the charging 
energy of the electron ground state). At higher gate voltage, Ee + aVa — f3Vg + Ue moves below 
the emitter chemical potential, and A^e approaches 1. When electron- hole pairs are generated due 
to optical pumping in addition to the tunnelling electrons, various exciton complexes such as the 
exciton, positive and negative trions, and biexciton can be formed. These new channels allow carriers 
in the emitter electrode to tunnel into the QD at lower gate voltage. When the gate voltage is below 
Vg2, but above Vgi, only the positive trion state is below the emitter chemical potential, hence only a 
small plateau is seen. As the gate voltage increases, more exciton complex levels are depressed below 
the emitter chemical potential, and a series of plateaus appear. The value of A^e depends on which 
exciton complex state is filled and on the pumping intensity. 

The tunnelling current as a function of the gate voltage can be measured directly. We see that 
optical pumping leads to two additional peaks below the threshold voltage Vg^, which is caused by 
the electron tunnelling assisted by the presence of a hole in the QD. This interesting phenomenon 
was observed by Fujiwara et al. in an SET composed of one silicon QD and three electrodes. The 
behavior of the photo-induced tunnelling current can be understood by analyzing the energy poles 
of the retarded Green's function as given in Eq. (12). In Figure 4, the first peak of the dashed hne 
is caused by a resonant tunnelling through the energy pole at e — E^ — 2Ueh (which corresponds 
to a positive trion state) with probability (1 — Ne-a)i'>^h,<Tnh-a)- The second peak is caused by a 
pair of poles at e = Ee — Ueh (the exciton state) and e = Ee + Ue — 2Ueh (the biexciton state) with 
probabilities (1 - Ne-a){nh,a + Uh-a - 2nh,anh-a) and Ne-a{nh,anh-a)- Since the magnitude of Ue 
is very close to that of Ueh for the present case, the corresponding two peaks merge into one. For 
other cases (e.g. smaller QDs), the difference between Ue and Ueh may be large enough for the peaks 
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to be distinguished. The third peak is caused by another pair of poles ed, e — (the single-electron 
state) and e — Eg + Ue — Ueh (the negative trion state) with probabilities (1 — Ne,a){i — {'nh,a + 
rih-a) + nh^anh-a) and N^-ainh^^ + rih-a — '^nh^a'^h-a)- The last peak locating near Vg = 123mV 
is due to the resonant tunnelling through the energy level at the pole e = E^ + Ue (the two-electron 
state) with probability iVe,_o-(l — (n/i^CT + ^^/i.-fj) +'T'/i,o-'^/i,-o-)- Obviously, the relative strengths of these 
tunnelling current peaks are influenced by the hole occupation numbers. The gate voltage differences 
^^21 = Vg2 — Vgi ^ud V^43 = Vg4 — Vg3 allow the determination of the electron-hole interaction and 
the electron-electron Coulomb repulsion, since /3AV^2i = Ueh and /3lp43 = Ue- Consequently, such a 
measurement can be used to determine the exciton binding energy, which is otherwise not possible 
via typical photoluminescence or the electrically driven emission spectrum measurement. In such a 
measurement, it is necessary to distinguish Ueh from Ue- Thus, we recommend using a smaller bias 
(V^), which would make the tunnelling peaks sharper and the double-peak features at = Vg2 and 
Vg — Vg3 better resolved. 

Figure 5 shows the tunnelling current as a function of gate voltage for various strengths 
of photo-excitation power at zero temperature. According to the definition of ■ye{h),c = 
- Eq \geih),ci\'^Im l/{-fe(h),s{% + heih),s)), the coudition = T/j > 7e,c is still satisfied, even though 
we tune / = 'yh,cNh,k/Th,s up to 10 in Fig. 5. 7=10 can be regarded as increasing pexc by 10 times, 
but it is still in the weak pumping regime. We notice that increasing photo-excitation power tends to 
suppress the tunneUing current at high gate-voltage. For instance, the peak located at Vg — 123mV 
almost vanishes due to the much reduced probabihty for the pole Bit e — Ee + Ug. Based on the results 
of Figs. 4 and 5, we see that the SET has potential apphcation as an optically controlled switch. For 
practical applications, we must consider the behavior of the tunnelling current for an SET operated 
at finite temperatures. Figure 6 shows a comparison of the electron occupation number and tun- 
nelling current at zero temperature and finite temperature (/csT — 2meV). We see that the plateaus 
become broadened and the magnitude of the tunnelling current peak is reduced as the temperature 
increases. Thus, a large charging energy is required for SETs to be operated at high temperatures. 
The tunnelling current as a function of applied bias (14) for various strengths of photo-excitation 
power at ksT = 2meV and Vg = is shown in Fig. 7. The tunnelling current exhibits a staircase 
behavior. This is the well-known Coulomb blockade effect. The effect of the optical pumping is reduce 
the threshold voltage and increase the number of plateaus in the J-V characteristics. We also notice 
that a negative differential conductance occurs at high voltages when the resonance level of the QD 
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is below the conduction band minimum of the emitter electrode. 

B. EfTects of electrodes on the spontaneous emission spectrum 

SET has been considered as a single-photon generater (SPG) for quantum information application 
such as quantum cryptography and teleportation.^^ Thus, it is of great interest to study the sponta- 
neous single-photon emission spectrum of the SET. We previously reported the theoretical studies of 
the emission spectrum of a single QD embedded in a p-n junction. ^^'^^ Here, we consider the single- 
photon emission of an SET (a QD embedded in a n-i-n junction) under optical pumping. Since the 
magnitudes of inter-particle Coulomb interactions Ug, Uh and Ueh are fairly close, it will be difficult 
to resolve the peaks of exciton ( X), negative trion (X^), positive trion {X~^) and biexciton (X^) for 
the case of large tunnelling rate as we considered above (F^ = = 0.5meV). Therefore, for this 
study we adopt a small tunnelling rate, Tl = Tr = O.lmeV in order to resolve these peaks in the 
emission spectrum. This can be achieved, for example, by increasing the barrier thickness between 
the QD and the conducting leads. Figure 8 shows the emission spectrum for various gate-voltages 
at fixed pump power (/ = 0.9). For Vg — lOOmV, the calculated spontaneous emission spectrum 
displays four peaks corresponding to X~, X, X"^ and The red shift of the X~ peak relative to 
the exciton peak agrees well with experimental observations. The biexciton peak displaying 
a blue shift with respect to the exciton peak (showing an antibinding biexciton) is also consistent 
with the observation reported in ref.[14, 19-21] Recently, studies of the binding and antibinding of 
biexcitons were reported in Ref. 21. Our calculation given by Eqs. (21) and (22) provides only the 
antibinding feature of biexciton. In Ref. 21 it is pointed out that the biexciton complex changes from 
antibinding to binding as the QD size increases. For QDs with dimension larger than the exciton Bohr 
radius (around 20nm), the correlation energy becomes significant. In this study, we have not taken 
into account the correlation energy. This is justified as long as we restrict ourselves to QDs with size 
less than 20nm. We see in this figure when the gate-voltage is increased up to = 12bmV, only the 
X~ and X"^ peaks survive in the emission spectrum as approaches 1. Furthermore, we can adjust 
the intensity of excitation power to either enhance or suppress the strength of the biexciton peak. 
Thompson et al. demonstrated that the single photon generated by the biexciton state is significantly 
less in emission time than the single exciton state. From this point of view, the biexciton state is fa- 
vored for the application of single-photon generation. We find that one can increase the pump power 
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in order to select the biexciton state as the desired source to produce the single-photon emission. 
This is demonstrated in Fig. 9, where the emission spectra for various strengths of excitation power 
at the fixed gate voltage Vg = 125mV are shown. We see that the negative trion saturates at / = 1.9 
and diminishes afterwards, while the biexciton peak increases quadratically with / and becomes the 
dominant peak in the emission spectrum. 

Finally, we discuss the emission spectrum of an isolated charge- neutral QD with = = 0. 
For this case the rate equation of exciton number Nx,a can be readily obtained by solving the rate 
equation 

7x,ciVfc(l - Nx,a) f —up{u)\ImP{u)\ = (24) 

J TT 

with 

where f2 = 47rn^/i^/(6c^^^eo) and is the refractive index of the system. p(uj) = u'^ arises from the 
density of states of photons. Ex and AEx,c denote, respectively, the exciton energy and exciton- 
exciton correlation energy. 7x,c denotes the capture rate for excitons from the wetting layer to the 
QD. After the integration of u; in Eq. (24), we obtain 

QNxAC^ - Nx,-a)E^x + Nx,-a(Ex + AE,f) (26) 

for small F, and 

where 7x,r = ^Ex denotes the exciton decaying rate, and we have ignored AEc, since it is small 
compared to Ex- According to Eq. (25), we see that the intensities of exciton and biexciton peaks 
are proportional to Pexc and pl^^, if Nx,a is proportional to Pexc (the excitation power). This is in 
very good agreement with many experimental reports. ^^'^^"^^'^^ Note that in some experiments, the 
pump energy is far above the QD band gap and the emission occurs after the electrons and holes 
are captured. This situation becomes similar to the open system considered in the SET with optical 
pumping. In this case, we can replace N^^a in Eq- (20) by rie^a, the electron occupation due to optical 
pumping, and we can still prove that the intensities of X and peaks as functions of the excitation 
power are described by Eqs. (14), (15) and (20). Thus, they exhibit linear and quadratic behavior as 
a function of Pexc, respectively. So, there is no qualitative difference between Eq. (20) and Eq. (24) 
regarding the behavior of X and X^, except that Eq. (20) also describes the possibihty of emission 
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due to recombination in trions. This provides an explanation to the experimental observation of Ref. 
17, in which the pump energy is far above the QD band gap and both X~ and peaks (in addition 
to the X and X'^ peaks) were observed in the emission spectrum. 

VI. SUMMARY 

In this article we have studied the tunnelling current and emission spectrum of an SET under optical 
pumping theoretically. We apply our theoretical analysis to an SET, which consists of a single quantum 
dot embedded in a n-i-n junction. We use a simple but realistic effective mass model to calculate the 
inter-particle Coulomb interactions of pyramid-shaped InAs/GaAs QDs. These Coulomb interactions 
are essential for the study of optical and transport properties of QDs. The retarded Green's function 
was calculated via a non-equilibrium Green's function method. It is found that the holes in QD 
generated by optical pumping leads to new channels for electrons to tunnel from the emitter to the 
collector, which creates a variety of ways to control both the electrical and optical signals. The binding 
energy of exciton complexes as well as electron charging energy can be determined by examining the 
tunnelling current as a function of gate voltage Vg with a small applied bias V^. In addition, the 
emission spectrum of SET can be modified significantly by adjusting bias voltage, the gate voltage, 
or the pump intensity. 
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VII. APPENDIX A 

In the system described by Eq. (1), electromagnetic field of frequency cuq is used to resonantly 
pump carriers into the InAs wetting layer. Therefore, electrons can be injected into the InAs QD 
via either the tunnelling process or carrier capture from the wetting layer. On the other hand, holes 
in InAs QD are provided only by the capture process. The dominant carrier capture process is the 
the optical-phonon assisted carrier relaxation as described by the fifth term in Eq. (1). Because 
of the phonon-bottleneck effect, the carrier capture rates for QDs are suppressed. Consequently, 
the condition Fg ^ 7e,c is usually satisfied, where and 7e,c denote the electron tunnelling rate 
and electron capture rate, respectively. Based on the above condition, we can ignore the electron 
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capture process in the calculation of the tunnelling current. The tunnelling current from the left 
electrode to the InAs QD can be calculated from the time evolution of the occupation number, 

Before the calculation of tunnelling current, an unitary transformation is employed to remove the 
phase of the optical pumping term. The renormalized energy levels of the electron and hole become 
Ee — Ee — UJ /2 and eh — — uj /2. Furthermore, the hopping terms between the leads and the dot 
become time dependent, Vkj(i) = Vkexp[—{-^iujt) /2]. Now, we can apply the formulas given in ref.[5]. 
We obtain 

J=^j ^rL(6)[iG<,(6) + fL{e)ImGlSe)l (28) 

where Ti,{e) — 27rX^k |yL,kp(^(e — ^k), and fhi^) is the Fermi distribution function. Gf^^{e) and Gl^^{e) 
are the Fourier transformation of the lesser Green function Gf^^{t) — i(d^^(0)de,o-(0) and retarded 
Green's function Gl^^{t) = —i9{t){{de^a^{t),dl^^{0)}). Applying Dyson's equation, we obtain 

GU')^r{e)[GlAe)-GlM (29) 
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with 



where gf^^^^.^{e) = h{e){gl^^^.^{e) - <k,a;^(e))-(i/I(e) - 9i{^)) = «27r(5(e - et). The lesser Green's 
function and retarded Green's function for holes are denoted by Q^-d^) ^'^^ gti-ai^)- 
small electron-hole recombination rate R^h ~ terms involving |Ap can be droped. That is to 

ignore the photocurrent in Eq. (28). After rewriting /^(e) as 

^ TL{e) + Vn{e) ' ^^^^ 

and substituting it into Eq. (28), we obtain the tunnelling current 

The expression of Eq. (32) was also obtained by Meir, Wingreen and Lcc^. It is of no surprise 
that we reproduce their result, because the system considered here reduces to their case as the carrier 
relaxation process and electron-hole recombination process are both turned off. Note that the retarded 
Green's function Gg ^.(e) includes the full effects due to the electron-electron Coulomb interaction and 
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electron-hole Coulomb interaction and the effect due to the dot-electrodes coupling to the lowest 
order. Also note that the electrodes are coupled to the conduction band of the QD, but not to the 
valence band of the QD. 

VIII. APPENDIX B 

When light is used to resonantly create electrons and holes in the wetting layer, electrons and 
holes can relax to the InAs QD via optical phonon-assisted process. When electrons relax into InAs 
QD, they will tunnel out quickly since the electrons in the ground state of the InAs QD are coupled 
with the electrodes. However, holes in the InAs QD can decay only via recombination with electrons 
(radiation process) or via impurity scattering (nonradiation process). When holes are present in 
the InAs QD, we expect the tunnelling current to change significantly due to the large electron-hole 
Coulomb interactions [see Eq. (12)]. In this appendix, we give the derivations for Eqs. (14) and (15), 
which describe the electron and hole populations in the InAs QD under optical pumping. First, we 
introduce the electron and hole distribution functions in the wetting layer, A^g k = (^e,k(0^e,k(^)) and 
Nh,-k — {dl _i^{t)dh-k{t)). We use the equation-of-motion method to calculate A^e,k = (^I,k(^)^e,k(^)) 
and Nh-k = (4 -k(^)4 -k(t))- We obtain 

^A^e,k = -hcsNeM + 2/m(Aj6^4,-k) 

-2^/m5e,q(^&UO, (33) 
q 

and 

^^^/^,-k = -^7/.,s^^/^,-k + 2/m(Aj6lX-k) 

-2j2lmg,,^{Abl_^dn,-a), (34) 
q 

where 7e(h),s denotes the carrier scattering rate due to mechanisms not considered in the Hamilto- 
nian. (Aofee,k^I,-k) represents the polarization for creating an electron-hole pair in the wetting layer. 
ge,ci{Ablkde,a) (fl'/i,q(^^I,-k^/i,-o-)) tcrms due to the process for transferring one electron (hole) 
from the QD to the wetting layer by absorbing one optical phonon. To terminate the equations of 
motion [Eqs. (33) and (34)], we need to solve (Ao6gk^I,-k) ge,q{Abli^de,a) (5'/i,q(^^I,-k'^/i,-(j))- 
The equation of motion for ge,q{Abl^^de,a) leads to 

^9e,ci{Abl^de,a) = {i^q - 7e,.) (^&I,k4,a) + IQcmI^ 
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[(1 + n,)(l - K + A^e))A^e,k - + ^e)(l - A^e,k)], (35) 

where denotes the phonon distribution function, defined as nq{u!q) = l/{exp{ujq/ki)T) — 1). He = 
{dl^^{t)de,(^{t))p and A^e = {dl^a{'t)de,a{'t))t denote the electron occupation number of the QD due to 
the pumping process and tunneUing process, respectively, fig — (ek — Eg — ujq). Meanwhile, the 
polarization obeys the following equation 

-i|Aor(l-iVe,k-A^/.,-k), (36) 
where flo = ee,k + e/i,-k — i-^o- In the steady state, we obtain 

Ne,k = Pexcil - Ne,k - Nh,-k) - 7e,c[(l + nq){l - (Ue + A^e)) 

A^e,k-n,(ne + iVe)(l-iVe,k)], (37) 

where Pexc = -Im\Xo\^/ {^e,s{^o + i{le,s + lh,s))) and 7e,c = - Eq Im\ge,ci\'^ / {'^e,s{^q + ^le,s))- For low 
pumping intensity, A^e,k is linearly proportional to Pexc- In addition, we note that 7e,c is fairly small 
due to the phonon bottleneck effect of the QD. Similarly for holes in the wetting layer, we obtain 

Nh-ii = Pexc(l - A^e,k - A^/i -k) - 7fe,c[(l + ^9)(1 " ^h) 

Nh,-^.-nqnh{l-Nh,-^,)l (38) 

where ^h,c — —H(ili^\gh,(i\^ /{lh,s{^q + ilh,s)) is the hole capture rate. Next, we need to solve the 
equation of motion for ne,c and Uh-a- We have 

+ 2/m<7e,q(A6U4,a), (39) 

and 

+ 2/m^ft,q(^4^_k4,-a), (40) 

where — -\-Vr and F^g denote the electron tunnelling rate and hole- impurity scattering rate, 
respectively. At low temperatures, the phonon-absorption process is negligible, since Uq — > 0. Using 
Eqs. (35) and (20), we obtain the steady-state solution for ne,a 

ne,. = 7 V, 41 

7e,cJVe,k + RehUh + 
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Note that the electron-hole recombination rate R^h is given by 

Reh^^l —p{uj)\ImV{u;)\, (42) 

J TT 

where Q = 47rn^/i^/(6c^^^eo) and Ur is the refractive index of the system. p{uj) — u? arises from the 
density of states of photons. Similarly, we obtain 

Because of the condition ^ 7e c and weak pumping intensity, we have A^j^k ~ Pexci and Ue is 
neghgible. As for N^, we can employ Eqs. (29) and (31) to obtain 

de TLfL{e) + TRfR{e) 



.--f- ^'^fll^'^^^'^ IrnGlAey (44) 
J n i L + i R 
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Figure Captions 

Fig. 1: Schematic diagram of a single electron transistor under optical pumping. Fl and Fr 
denote, respectively, the tunnelling rate for electrons in the QD to the source and drain electrodes. 
Reh denotes the electro-hole recombination rate. The vertical wavy line indicates the optical pumping, 
and the horizontal wavy line indicates the single-photon emission. 

Fig. 2: The lowest three energy levels of a quantum dot as functions of the QD size b for (a) 
electrons and (b) holes. 

Fig. 3: Intralevel Coulomb interactions C/g and Uh and inter level Coulomb interaction Ueh as 
functions of the QD base length b. 

Fig. 4: Electron occupation number N^. and tunnelling current as functions of gate voltage at zero 
temperature for various strengths of optical excitation. Current density is in units of Jo = 2ex meV/ h. 

Fig. 5: Tunnelling current as a function of gate voltage at zero temperature for various strengths 
of optical excitation. 

Fig. 6: Electron occupation number Ng and tunnelling current as functions of gate voltage for 
various temperatures at fixed optical excitation strength (/ = 0.9). Current density is in units of 

Jo = 2e X meV/ h. 

Fig. 7: Tunnelling current as a function of applied bias at temperature ksT = 2meV for various 
strengths of optical excitation. Current density is in units of Jo = 2e x meV/ h. 

Fig. 8: Intensities of emission spectrum for various gate voltages at fixed optical excitation strength 
(7 = 0.9). 

Fig. 9: Intensities of emission spectrum for various strengths of optical excitation at fixed gate 
voltage {Vg = 125mV). 
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